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ABSTRACT 


The problem considered was the development of methods to 
mathematically model major terrain features for an autonomous 
surface vehicle guidance system. The terrain information 
available was restricted to data that could be measured by the 
vehicle's terrain sensor, which was assumed to provide range 
readings for specified values of azimuth and elevation angles. 

Terrain indices were developed to process the sensor 
measurements. The purpose of these indices is to enable the 
vehicle to directly or indirectly detect major terrain features 
that represent an obstacle to travel. Modeling methods were 
then constructed from the indices to provide a systematic means 
of modeling the terrain. 

Simulated sensor measurements were obtained by using the 
geodetic survey map of a relatively hilly portion of terrain. 
The simulated sensor data was then used to evaluate the effec- 
tiveness of the modeling methods. The effectiveness of a 
particular method was- primarily defined in terms of the per- 
centage of passable terrain that could be detected and the 
error in range due to overestimating the distance to a terrain 
obstacle. 

It is concluded that it is feasible to model major terrain 
features for a vehicle guidance system with only a modest on- 
board computer requirement. Three alternative modeling methods 
were developed that are applicable to relatively rough terrain. 
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The method chosen for a particular situation would depend on 
the modeling effectiveness desired and the sensor and computer 
capability available. 
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I . INTRODUCTION 

The problem under consideration is the mathematical 
modeling of major terrain features using line-of-sight sensor 
measurements acquired from a single terrain location. This 
investigation is part of an overall National Aeronautics and 
Space Administration sponsored study of a surface vehicle con- 
trol system for the unmanned exploration of Mars, It is 
assumed that the exploration vehicle will have to travel a 
considerable distance, perhaps a hundred kilometers or more, 
from its initial landing site to one or more remote locations 
to complete its mission. The vehicle is considered to be oper- 
ating independent of Earth-based control or other orbiting 
satellites. Thus, the guidance system must be self-contained, 
arid it must be capable of processing terrain information and 
making a path selection decision. 

The path selected should not only lead the vehicle toward 
its objective, but also should provide a safe path of travel for 
the vehicle. If the vehicle should encounter a terrain hazard, 
the mission could end in failure. Thus, the guidance system 
must be able to: 1). avoid short range obstacles that represent 
an immediate danger to the safety of the vehicle, and 2). select 
a long range path around major terrain features that obstruct a 
direct and safe route to the objective. This report deals 
specifically with the terrain model required to perform long 
range path selection. 
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A number of path selection algorithms (Ref. 1,2) have 
been developed which can be used for long range vehicle guidance. 
In order to utilize these algorithms, however, the terrain must 
be modeled with respect to the obstacles that are in the vicinity 
of the vehicle location. The terrain model required is two 
dimensional in nature and simply describes the distance a 
vehicle can travel before encountering an obstacle in a series 
of discrete azimuth directions. A sample terrain model is 
pictorial represented in Figure 1, 

Basically, there are three main steps that must be performed 
by a terrain modeling system, namely: 

a) Terrain Data Acquisition 

* 

b) Terrain Data Processing 

c) Terrain Model Generation 

Section II deals specifically with terrain data acquisition. 
Since specifications of a sensor to acquire terrain data have 
not been defined, typical long range sensor measurements are 
simulated to indirectly determine some general sensor measure- 
ment characteristics. The sensor data simulation is performed 
by assuming various vehicle locations on a geodetic map, and 
then converting topographic information into sensor data, i.e., 
range measurements for specified values of azimuth and ele- 
vation angles . 

In Section III, a number of indices are developed to pro- 
cess the sensor terrain data in order to detect the presence 
of a terrain obstacle. The terrain indices are combined in 




FIGURE 1. Terrain Model 




Section IV to provide several systematic methods to generate a 
terrain model. The simulated sensor data is then used to 
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evaluate the effectivenss of the modeling methods. The effective 
ness of a method is primarily defined in terms of the per- 
centage of passable terrain that can be detected and the error 
in range due to overestimating the distance to a terrain 
obstacle. Then, some general sensor measurement requirements 
are described in terms of the effectiveness of the modeling 
methods. 

II. TERRAIN SIMULATION 
A. Sample Terrain 

There are four general terrain characteristics that 
significantly effect a vehicle Vs performance (Ref. 3,4,5). 
These characteristics are: 

a) surface composition 

b) surface geometry 

c) vegatation structure 

d) hydrologic geometry 

Presently, .there is relatively little known about 
the surface conditions on Mars. From the general infor- 
mation that is available, there is no evidence of vega- 
tation or significant amounts of water present. The sur- 
face composition is a matter of speculation. Considering 
the surface geometry, it is estimated (Ref. 6) that slopes 
of major terrain features generally vary between 0 and 30 
degrees, and that over half of the slopes are less than 
10 degrees. 
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Thus, the surface geometry is the only terrain character- 
istic which will be considered in this study as con- 
stituting an obstacle to travel. 

For the purposes of long range terrain modeling, a 
slope of a major terrain feature equal to or greater than 
10 degrees will be defined as an obstacle. Slope, in its 
classical sense, is defined as the angular deviation of a 
surface from the horizontal, measured perpendicular to 
the topographic contours. Slope so defined can be 
measured from a contour map. A major terrain feature is 
defined as an area that is significantly larger than the 
dimensions of the exploratory vehicle (Ref. 7), The 10 
degree slope, as the definition of an obstacle, was some- 
what arbitrarily chosen, and it may even be considered 
conservative with respect to the climbing capabilities of 
proposed exploratory vehicles. In justification of this 
obstacle definition, it should be kept in mind that an 
exploratory vehicle may encounter adverse soil conditions 
or some unknown element that could substantially reduce 
the vehicle's climbing capabilities. 

Slope was the only surface geometry factor con- 
sidered as a true obstacle. Such factors (Ref. 4,8,9,10) 
as slope segments, spacing and size of vertical obstacles, 
terrain approach angles, etc, are pertinent to short 
range terrain modeling and will not be discussed further. 
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A geodetic map (Ref. 11) was used to represent a 
sample terrain and was the source of empirical terrain 
data for the terrain modeling techniques developed. The 
map selected represented a section of hilly terrain near 
the Catskill region in New York State with a scale of 
1:24000 and a contour interval of 10 feet. The slopes of 
major terrain features varied between 0 and 30 degrees, 
whereas the slopes of typically encountered terrain 
features varied between 0 and 15 degrees. This slope 
variation is somewhat comparable to the values for Mars 
referred to earlier. 

Several significantly different geometrical terrain 
configurations were specifically analyzed in detail in 
the hope that the results of this investigation would be 
applicable to a variety of terrain situations. The 
three areas selected are listed below with a brief descrip- 
tion of their salient features, as seen from the 
assumed positions of the vehicle. The sectors are 
indicated in Figures 2 and 3. 

Area 1: The vehicle's vantage point is on the summit 

of a hill (reference map grid coordinate 16420962*). The 
vehicle is sighting east across a valley toward a ridge 
line. A 35 degree sector of terrain was considered (see 

* The first four digits refer to the east-west grid lines, and 
the last four digits refer to the north-south grid lines. 
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Figure 2) . 

Area 2: The vehicle's vantage point is on a gentle 

rise (reference map grid coordinate 14601673). The 
vehicle is sighting north-northwest toward relatively 
minor hills and valleys scattered through the area. A 50 
degree sector of terrain was considered (see Figure 3) . 

Area 3: The vehicle's vantage point is on a gentle 

rise (reference map grid coordinate 14601673). The vehicle 
is sighting north-northwest toward a relatively high moun- 
tainous area. A 50 degree sector of terrain was considered 
(see Figure 3) . 

In addition to the areas described above, other 
selected profiles of terrain were also considered, such as 
from a valley location sighting across a flat area towards 
a mountainous region, and from the summit of a hill sight- 
ing toward a valley. 

B. Terrain Data Acquisition 

The acquisition of data by a terrain sensor on the 
vehicle was simulated by making measurements on the geo- 
detic map in the areas decribed above. 

Three types of sensor measurements, i.e., azimuth 
angle, elevation angle and range, are assumed to be avail- 
able, from which all terrain information must be obtained. 
The azimuth angle ( ©) represents the angle the terrain 
sensor is directed in the horizontal plane of the vehicle. 
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Angles of Q will be measured with reference to some con- 
venient azimuth direction, e.g. , grid north on the geodetic 
map is 0°. The elevation angle ( (3 ) is the angle the 
sensor is directed in the vertical plane of the vehicle. 
Angles of will be measured with reference to the 
horizontal plane at the vehicle's location, with angles 
above the horizontal considered positive. The range (R) 
is the distance the sensor can measure along the line of 
sight to a terrain feature. The angles © and f3 are 
adjusted at the vehicle and thus are independent variables, 
whereas R is a dependent variable. It is assumed that 
values of R can only be measured at discrete increments 
of © and , and that the size of these discrete 

increments is a characteristic of the sensor which must 
be specified. 

In order to simulate the terrain data which the 
vehicle would receive from a particular vantage point, 
values of distance and relative altitude (Z) were measured 
directly from the geodetic map for a fixed value of 
azimuth angle (©). From this information, terrain pro- 
files of altitude vs. distance plots were constructed 
for each value of © as illustrated In Figure 4. 

The procedure used to obtain simulated sensor data 
was to assume a sensor range measurement, R, to a point 
on the terrain profile. Then the corresponding elevation 





angle, , 


could be calculated from 
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sin ^ - 


Z 

R 


( 2 . 1 ) 


Figure 5 is the corresponding f3 — R plot for the terrain 
profile presented in Figure 4. The points along the con- 
tinuous R curve are the discrete sensor data the vehicle 
would receive in the absence of measurement errors. It 
is noted that if areas of the terrain profile were hidden 
from the vehicle’s view, there would be a corresponding 
discontinuous horizontal segment on the R plot, as 
illustrated in Figures 4 and 5. 


Sensor measurements were also compiled in the form of 
® vs. R polar plots, with ^ held constant. The polar 
plot can be constructed from either the ^3 — R plots taken 
in a series of azimuth directions or directly from the 
geodetic map using templates graduated in Z and R for 
fixed values of ^3 • The O — R polar plots are useful in 
detecting major terrain formations but are only indirectly 
related to the slopes of these terrain formations. 

Figure 6 is a 0-R polar plot acquired with {3 fixed at a 
small positive value. The areas labeled with a C were 
relatively clear of high terrain formations, i.e. , no 
terrain feature was sensed in the ranges considered. Con- 
versely, the areas labeled A,B,D and E corresponded to 
relatively high terrain formations. 
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C. Characteristics of Sensor Measurements 

Sensor measurements of the sample terrain were pre- 
dominantly compiled in the form of ^3— -R plots. These 
plots were found to be more useful than the ^ — R plots for 
analyzing the sensor data and developing terrain modeling 
techniques since they yielded information directly 
related to the slope of the terrain along a profile. The 
( 3 — R plots were prepared in azimuth increments (AG) 
of 5 degrees for the areas of the map described above. A 
total of 35 such plots were drawn. Since the £3 — R plots 
were acquired from various types of terrain configurations, 
certain sensor measurement characteristics were deduced 
from the data. 

A range of 400 feet from the vehicle was the minimum 
value for which useful empirical data could be obtained 
from the geodetic map due to the lack of sufficient con- 
tour information in the proximity of the assumed vehicle 
locations. The vehicle locations were specifically 
chosen so as not to be in areas of large contour 
variation since. the purpose of the guidance system is to 
keep the vehicle away from such areas. Due to the lack 
of short range contour information, sensor measurements 
often had to be approximated when R was less than 400 
feet. Thus, any conclusions drawn from the simulated 
terrain data only pertain to terrain at ranges greater 
than 400 feet from the vehicle. 
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The maximum feasible range for modeling the sample 
terrain was chosen as 4000 feet since there were insufficient 
sensor measurements available to the vehicle at greater 
range values. Figure 7 is a h 1st or gram of the number of 
sensor measurements acquired from 35 R plots (2y 3= ,005 
radian) as a function of the range at which the measure- 
ments were sensed. Nearly 90% of the discrete sensor 
measurements were sensed at ranges of less then 3500 feet. 

The maximum value of needed to acquire sensor 
information was .073 radian (4.2 degrees). Approximately 


80% of the sensor measurements were acquired for values of 


^5 = +.030 radian (1.8 degrees). 


For all values of ^ 


encountered, small angle approximations were certainly 
valid . A typical use of this approximation was in the cal- 
culation of ( 2 > from the terrain profile measurements using 



( 2 . 2 ) 


III. TERRAIN DATA PROCESSING 


The terrain data in the form of a set of sensor measure- 
ments R^,^ and 0 ^ must be processed to yield information on 
the presence or absence of a terrain obstacle (slope). Each 
terrain index that is presented processes the sensor measure- 
ments in a particular manner to yield a numerical quantity 
which is compared to an index test number. The test of the 
numerical quantity is made to detect the presence of a terrain 



Number of Occurences 


L6 



R (feet) 

FIGURE 7. Histogram of the Distribution of Sensor Information 
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obstacle. 


All of the terrain indices use ^ vs. R sensor information 
that can be acquired along fixed azimuth directions (©j ). It 
is assumed that the sensor can acquire successive range measure- 
ments (R) along an azimuth direction for increasing values of 
elevation angle (^3 ). Thus, R will be a monotone increasing 
function along a terrain profile in all cases of interest. 

A. Slope Index 

A calculation of the linear approximation to the slope 
along a terrain profile is made. The slope along a pro- 
file is defined as the in-path slope (S^- ) and it is 
derived in the Appendix in terms of sensor measurements. 

The formula for Sj is 


S T = /02 R 2-(#1 R 1 (3.1) 

r 2 * R 1 

where the quantities^ R^ and^^j &2 denote successive 
measurements of elevation angle and range for a fixed 
azimuth. Generally ,^2 w ill be equal to ^3 ^ 4- A^> . The 
numbers "l" and "2" denote the first and second measure- 
ments of any pair of constant azimuth readings. 

The slope index compares the calculated value of Sj 
against the index test number, T s , by ascertaining if the 
inequality 


P Z S. -A R 

r 2 ' R 1 


1 T. 


(3.2) 



is satisfied. If it is, then is taken as the estimated 
range to an obstacle along that particular azimuth, © j . 

In the testing done, T s was taken as tan 10°. 

As the size of used to acquire terrain infor- 

mation is increased, the calculated value of Sj tends to 
yield a lower estimate of the true slope. In this case 
the linear approximation of the slope is necessarily 
made over a longer segment of the terrain profile. 

Thus, if a large value of is used, it may be necessary 
to use a lower value of T g in order to avoid local slopes 
which would exceed 10 degrees. 

The index would be of little merit if corresponded 
to a different terrain feature than R£ . In this case, 
the calculated value of S-j- would then only represent the 
slope of an imaginary line connecting the two terrain 
features sensed. A means of overcoming this situation 
is to use another terrain index in conjunction with the 
slope index that would be applicable for this situation. 

Gradient Index 

A numerical approximation to the gradient of a terrain 
feature is made by varying both © and ^ for a set of 
measurements. The formula for the gradient calculation 
is developed from the magnitude of the gradient (Ref. 12) 
in polar coordinates where the altitude Z is described as 



19 


a function of position, Z = Z(D,©), to yield 

“ I -JftsH* !§)' 


(3.3) 


The first term in (3.3) can be approximated by as 

indicated in the Appendix, calculated along the line of 

constant azimuth © = ©. . The second term in (3.3) can 

J 

be approximated by the cross-path slope (S c ) which is 
derived in the Appendix in terms of three sensor measure- 
ments acquired along a pair of adjacent azimuth directions, 
® j and The cross-path slope (S^ ) is the slope of 

the terrain feature that is perpendicular to the constant 
azimuth profile. The square root term in (3.3) may be 


expanded to yield 


G |= S T 1 + \ (f£) 2 - 1 (l£) 4 + 
1 2 S I 4 Si 


(3.4) 


provided that Sj is greater than S c . An upper bound 
on | G | is found by using only the first two terms of the 
expansion since the first term dropped is negative and 
dominates the error. 

The gradient index compares the calculated value of 
1 G against the index test number, T , by checking the 

■ O 

validity of the inequality 


S T + I > T 

i 2 Si § 


(3.5) 


If S G is found to be greater than Si , the square root 

(3.3) is expanded in terms of (S_ /S ) 2 . The 

JL C 


term in 



20 


gradient index would then involve the following comparison 


s ‘ + ^ T S 


(3.6) 


If either (3.5) or, if applicable, (3.6) is satisfied, then 
R i as given by (3.1), is taken as the estimated range to 
an obstacle along the line whose azimuth is ©^. 

The cross-path slope (S ) is calculated from sensor 
measurements made along the line © = © j and along the line 
corresponding to either © = O. + A© or © = ©.s~A©, If A$ 

J J 

is large with respect to the slope variation of the ter- 
rain, the value of S c calculated using©j and €? j + A© 


will differ substantially from the value calculated using 
& . and O. - hO . Thus, the effectiveness of the gradient 

V J 

calculation to detect the true slope of the terrain is not 
only a function of the size of A©, but also the choice 
of the adjacent azimuth direction used for the calculations. 


C. Relative Altitude Index 

A calculation of the relative altitude, Z, of a ter- 
rain feature is made in order to indirectly indicate the 
presence of a steep slope. The rationale for this cal- 
culation is that a steep slope is usually associated with 
a significant change in the relative altitude. The 
relative altitude can be calculated from one sensor 
measurement, using 

Z. = /S. R. 

l ( l l 


(3.7) 



The absolute value of the relative altitude is com- 
pared with an index test number, T r , by testing the 
following inequality: 




i R i 


> T, 


(3.8) 


If (3.8) is satisfied, an obstacle is assumed to begin 
at the previous range measurement along that azimuth 
line. The value of T r should be determined by the type 
of terrain that will be encountered, using either the 
statistical distribution of slope with Z or empirical 
test data. 

Since this test uses an indirect measure to find a 
steep slope it also should be stated that it could fail 
to identify a steep slope (i. e. , a steep slope not 
immediately corresponding to a significant change in Z) 
as an obstacle, or could label fairly safe but gently 
rising terrain as an obstacle. 

A significant feature of this index is the simple 
calculations required to detect a steep slope. Thus, 
this index should be relatively insensitive to sensor 
measurement ' errors . It is interesting to note that this 
index defines an obstacle in a similar manner as Lim 
(Ref. 2) for a myopic robot. 

Vertical Index 


A calculation of the change in successive values of 
R is made to indirectly indicate the presence of a steep 
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slope along a terrain profile. The vertical index com- 
pares the calculated change in R against an index test 
number, T v , and can be written as 

r 2 - \ £ T v (3.9) 

If (3.9) is satisfied, R^ is taken as the estimated range 
to an obstacle along the line © *> ^ . The rationale for 
this index is the following: since an increase in ^3 from 

to ^2 assoc: i- ate< 3 with a change in Z, a small change 
in R would correspond to a steep slope, as indicated 

in Figure 8. 

This index can be used in conjunction with the 
relative altitude index to partially overcome the limi- 
tation of the latter by detecting a steep slope that does 
not correspond to a significant change of Z. If T v is 
chosen to be a fixed value, it is very effective in 
detecting the presence of a 10 degree slope (for example) 
at one particular range, say R. However at ranges less 
than R the index is conservative since it considers 
slopes to be obstacles that are less than 10 degrees. 

Likewise, at ranges greater than R, it fails to detect 
a 10 degree slope as an obstacle. This effect is illus- 
trated in Table 1 where typical values of R£ - R^ are 
given if a 10 degree slope is to be measured (i.e.,2y3= .005 
radian) at various ranges from the vehicle. 
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TABLE 1 - Change in R for Measuring a 10° Slope 


R (feet) 

R 2 - R x (feet) 

500 

14.2 

1000 

28.4 

2000 

56.8 

3000 

85.2 

4000 

113.6 


If the formula for Sj is solved for the difference 
R 2 - R^, the resulting equation describes the variation 
of the range increments in terms of sensor parameters and 
Sj namely, 

R 2 ~ R 1 ^ 2 - R 2 (3.10) 

If T v were defined as the right side of (3.10), this 
index would become just a restatement of the slope index. 
Thus, T y will be chosen to be fixed. 

E. Discontinuity Index 

A calculation of the change in successive values of 
R is made to indicate the size of areas that are hidden 
from the vehicle for the obvious reason that these hidden 
areas may contain terrain obstacles. 

The discontinuity index compares the calculated 
change in R against an index test number, T^ , and can be 
written as 



R •> T, 

1 — d 


(3.11) 
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If (3.11) is satisfied, R ^ is taken as the estimated 

range to a potential obstacle along the azimith & =0 as 

J 

indicated in Figure 9. The value of T was obtained by 

d 

empirical data from the sample terrain. 

There are difficulties in assigning a value to T^ 
which is used to label hidden areas as obstacles. If T^ 
is made too large, there would be a high probability of 
encountering an obstacle in the hidden area. If T^ is 
made too small, relatively safe terrain would be labeled 
as being an obstacle. As an example of the latter, if a 
vehicle sensed a uniform sloping terrain of 1 degree, 
then, for A(^ - .005 radians, the difference R 2 - R^ 
would be only 100 feet at a nominal range of 1100 feet 
from the vehicle, but 1000 feet at a range of 3500 feet. 

IV. TERRAIN MODEL GENERATION 

A terrain modeling method is defined to be the repetitive 
application of certain terrain indices along successive azimuth 
directions to generate a table of azimuth values and ranges to 
the nearest obstacle. The following is a brief description of 
the operation of a modeling method: along each azimuth profile, 
the sensor measurements are tested until one of the indices 
indicates the probable presence of an obstacle. Then the 
range (R) to the obstacle and the corresponding value of O j 
are stored in a table in the computer memory. If no obstacle 
is found along a particular azimuth direction (©j), the last 
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R measurement and corresponding value of Q j are entered in 
this table. 

Alternate modeling methods are developed that differ in 
the type of computing capability required and the manner used 
to define terrain obstacles. The modeling methods formulated 
and tested are listed below, with their associated terrain 
indices : 


Method 1 : 


Discontinuity Index 
Vertical Index 
Relative Altitude Index 


Method 2: Discontinuity Index 

Slope Index 

Method 3: Discontinuity Index 

Gradient Index 

The first method detects obstacles by indirect means and 
requires little computing capability. The second method 
essentially uses the two dimensional terrain information along 
an azimuth profile to define an obstacle, while Method 3 uses 
three dimensional terrain information. Method 3 requires 
more computing capability than the other methods. 

The operation of the modeling methods was simulated using 
a WANG 370 Series Calculator (a programmable electronic desk 
calculator) . Methods 1 and 2 were programmed with 108 and 103 
operations (i.e.,add, multiply, change sign, etc.), respectively. 
Due to the limited program capacity. Method 3 was performed 
using two separate programs of 103 and 158 operations. 
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Storage and recall of previous sensor data was manually per- 
formed for this method. 

The simulation of the methods on an electronic desk cal- 
culator demonstrated the feasibility of using these same 
methods with a small onboard computer. The simulation also 
provided a practical means of processing large amounts of 
sample terrain data to test the effectiveness of the methods in 
modeling actual terrain. 

A. Comparison of Terrain Modeling Methods 

Terrain models of the selected areas described in 
Section IIA were constructed using each of the three 
modeling methods presented above. The modeling methods 
were designed to detect the passable terrain around the 
vehicle by processing sensor measurements. Passable ter- 
rain is defined in terms of the distance that a vehicle 
can travel in a given direction before encountering a 
terrain obstacle. In this study, the maximum possible 
range was limited to 4000 feet, and a slope in excess of 
10° was considered as an obstacle. Sensor measurements 
were acquired in increments of .005 radian (.3 degree) in 
elevation angle (&/&) and 5° degrees in azimuth angle (A©) 
to form a common basis to compare the individual modeling 
methods. The following index test numbers were used: T s = 
tan 10°, Tg = tan 10°, T r = 40 feet, T v = 25 feet, and 
T^ = 1000 feet. 



The effectiveness of a method to model terrain was 


evaluated by considering the following three factors for 
each of the three areas: 

1. The percentage of the passable terrain detected 
It is noted that terrain slopes of 10° or more 
contained within a 1° sector of terrain were 
considered as local terrain irregularities and 
not as terrain obstacles. 

2. The maximum range error computed. A range 
error is defined as the difference between the 
computed range (R) and the actual range to the 
first obstacle along a terrain profile. Only 
a computed R which is an overestimate of the 
range to the obstacle is considered. 

3. The percentage of terrain profiles along which 
range error was computed. 

The results of these tests are summarized below. 

Area 1: Method 1 found 39.77, of the passable ter- 

rain. The passable terrain that was not detected was 
the result of: 1). passable terrain was hidden from the 
vehicle's view, 2). some of the significant changes in Z 
did not correspond to a terrain obstacle. Due to the 

t 

convex shape of the terrain profile along two azimuth 
directions, large portions of passable terrain were 
hidden from the vehicle and thus were considered as 
obstacles. The existence of this hidden area also 
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effected the results obtained with the other two 
methods. A range error was found along 37% of the terrain 
profiles . The maximum range error was 150 feet which 
occurred for an actual obstacle at a range of 3250 feet. 
This range error was generally due to either a steep 
slope that did not correspond to a significant change in 
Z (i.e., T r ) or a change in R that was less than T v . 

.Methods 2 and 3 found 61.9% of the passable terrain. 

A range error was found along 25% of the terrain profiles . 
The maximum range error was 50 feet which occurred for an 
actual obstacle at a range of 3250 feet. The range errors 
were due to the value of Af 3 not being small enough to 
detect a narrow band (50 feet) of steep slope at a range 
of 3250 feet. 

Area 2: All three methods found only 9.2%, of the 

passable area. In this sector, relatively large areas of 
terrain near the assumed vehicle location were hidden 
from view. The discontinuity index was the primary means 
of detecting obstacles in this region.- No range errors 
were found in this sector. This was primarily because 
empirical data was used in selecting the value of T^. 

Area 3: There was a high percentage of passable 

terrain detected using all the methods since the assumed 
vehicle location offered a fairly good view of the ter- 
rain in the area. Method 1 found 61.5% of the passable 
terrain. The relative altitude index was effective in 
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detecting the presence of steep slopes because of the 
high correlation between the change in altitude and a 
steep slope. A range error was found along 27% of the 
terrain profiles. The maximum range error was 90 feet 
.which occurred for an actual obstacle at a range of 
1110 feet. The range error was due to the indirect 
means used to estimate the presence of a steep slope. 

•Method 2 found 80.4% of the passable terrain. A 
range error was found along 45% of the terrain profiles. 
The maximum range error was 390 feet which occurred for 
an actual obstacle at a range of 1000 feet. The large 
frequency and size of range errors was due to the fact 
that the slope along many of the terrain profiles was 
much less than the actual slope of the terrain, i.e., 
the cross -path slope was significant. 

Method 3 found only 54% of the passable terrain. The 
relatively low percentage of passable terrain detected was 
due to the "large" size of ^ & in relation to the size 
of the slope irregularities of the terrain. Both of the 
alternative azimuth directions were used in the gradient 
calculation. The most conservative estimate for the 
range to an obstacle was used for the terrain model. A 
range error was found along only 9% of the terrain pro- 
files, The maximum range error was 70 feet which occurred 
for an actual obstacle at a range of 2800 feet. An 



improvement in the range accuracy could be expected if a 
smaller increment in azimuth were used to acquire sensor 
measurements . 

Even though the effectiveness and applicability of 
each method varied with the terrain modeled, the follow- 
ing properties of the modeling methods were observed as 
a result of the above comparison. 

Method 1 detects a relatively small amount of pass- 
able terrain with modest size range errors. The effective 
ness of this method is primarily limited by the 
indirect means used to detect the presence of a steep 
slope. 

Method 2 detects the highest percentage of passable 
terrain but it results in relatively significant range 
errors if there is a significant cross-path slope 
present. 

Method 3 detects a modest percentage of the passable 
terrain with the size of A © used, but it results in the 
lowest range error of any method. The effectiveness of 
this method should be significantly increased if a 
smaller increment in azimuth angle was used. 

B. Terrain Sensor Measurement Requirements 

The effectiveness of the terrain modeling method 
depends on the characteristics of the terrain sensor 
available. A summary is given of some of the sensor 
measurement characteristics necessary to use the modeling 
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techniques developed, 


A terrain sensor was assumed to be available that 
could provide range measurements for distances up to 4000 
feet from the vehicle. The upper limit of 4000 feet was 
taken as the maximum range at which the modeling methods 
could be effectively applied to the sample terrain since 
sufficient terrain (sensor) information was not available 
at ranges greater than 4000 feet^as previously illustrated 
in Figure 7. There are two explanations why terrain 
information is more difficult to acquire as larger values 
of R are considered. First, terrain features relatively 
near the vehicle tend to obstruct line-of -sight measure- 
ments of terrain features relatively far from the vehicle. 
Second, the 1/R factor in the ^3 vs. Z relationship (2.2) 
tends to enlarge the minimum magnitude of Z that can be 
measured as R is increased, i.e. , for a fixed value of . 

If the sensor range measurements (R)are measured with 
some inaccuracy,- the effectiveness of a modeling method 
will be reduced corresponding to the type of calculations 
that are made with these range measurements. As an 
example, the slope index requires differencing, multi- 
plication and division operations, while the relative 
altitude index requires only the multiplication operation. 


The variations in the calculated values of the in-path 
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slope, e.g., S-j- = tan 10° - .176, and the relative 
altitude, e.g., Z = 40 feet for deviations of + 10 feet 
in the range measurements are shown in Table 2. 


TABLE 2. Effect of Range Measurement Errors of 
+ 10 Feet on Calculated Values of Z and 
S j (|3 l= 0 and 0 .01 radian) 


Nominal Range 
(feet) 

Calculated Z 

Jfeet} 

Calculated St 

Low 

High 

Low 

High 

1000 

39.6 

40.4- 

.130 

.276 

2000 

39.8 

40.2 

.152 

.224 

3000 

39.9 

40.1 

.156 

.200 

4000 

39.9 

40.1 

.162 

.194 


The maximum error in evaluating Z was 1% while the maxi- 
mum error in evaluating S-j- was 36%. Thus, Modeling 
Method 1 which includes the relative altitude index is 
less sensitive to range measurement errors than Method 2, 
which includes the slope index. It is noted that tan 10° 
and 40 feet were the values of T s and T v used in the 
simulation of the methods above to estimate the presence 
of an obstacle. 

2 . Elevation Angle Increment 

The size of the { 2 > increment used in acquiring sensor 
data directly affected the performance of the terrain 
modeling techniques. In Table 3, the percentage of the 
total passable terrain detected when the results for 
areas 1,2, and 3 are combined is tabulated as a function 
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of the increment AjO 


used to acquire sensor data. 


TABLE 3. Percentage of Passable Terrain 
Detected as a Function of A 


A ft 

(radians) 

Method 1 
(% passable 
terrain) 

Method 2 
(% passable 
terrain) 

.005 

31.7 

44.5 

.010 

32.1 

37,3 

.015 

17.8 

21.0 

.020 

9.8 

11.6 

.025 

8.2 

9.9 

.030 

7.3 

8.7 


There was a general decrease in the percentage of pass- 
able terrain detected when larger increments of ^<3 were 
used. This decrease in percentage was related to the 
amount of information that was available to the terrain 
modeling method. In Table 4, the average number of 
sensor measurements available was tabulated as a function 
of the size of i.e. , for a fixed azimuth direction. 

TABLE 4. Average Frequency of Sensor 
Measurements as a Function of &P 


Increment 

(radians) 

Average Number of 
Sensor Measurements 

.005 

8.2 

.010 

4.2 

.015 

3.0 

.020 

2.4 

.025 

2.0 

.030 

1.8 
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It would be possible to increase the percentage of 
passable terrain detected, if the index test numbers were 
optimized for the value of A/3 used. For example, the 
slope calculation of Sj yielded a lower approximation to 
the true slope if A^3 was large. Thus, a smaller value 
of T s should be chosen to account for the "cruder" slope 
approximation. 

3. Azimuth Angle Increment 

The size of the O increment used to acquire the 
sensor data should be related to the general character- 
istics of the terrain to be modeled. The increment in 0 
must be small enough to detect significant information 
about major terrain features that would be typically 
encountered . 

In general, the effectiveness of a modeling method 
decreased as the size of AO was increased. This trend 
was observed from the simulation work when the percentage 
of total passable terrain was seen to decrease as 
was increased, refer to Table 5. 

TABLE 5. Percentage of Passable Terrain Detected as a 
Function of AO = .005 radian) 


jA & 

Method 1 

Method 2 

(degrees) 

% passable terrain 

% passable terrain 

5 

31.7 

44.5 

10 

21.7 

26.7 
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V. DISCUSSION 

Three terrain modeling methods were developed in this 
report for an exploratory vehicle's long range path selection. 
The modeling method selected for a particular guidance system 
would depend upon the trade-offs made between the modeling 
effectiveness desired and the computer and sensor capabilities 
available. 

The performance of the methods was simulated on a sample 
terrain. The modeling methods could be adapted for terrain 
dissimilar to the sample terrain by proper choice of the ter- 
rain index test numbers. 

The main restriction on the terrain modeling effective- 
ness was the line-of-sight measurement requirement. If the 
effectiveness of a method was not acceptable, provisions might 
have to-be made to supply the vehicle with additional terrain 
information. Possibilities would include a second sensor 
located on a vertical extension of the vehicle which could be 
raised when the vehicle was stopped or perhaps data sent from 
a satellite orbiting overhead. 

There were several difficulties in interpreting discrete 
line-of-sight sensor measurements. For example, consider the 
two sensor measurements illustrated in the terrain profiles 
shown in Figures 10A and 10B. The same sensor measurements 
can be obtained from quite different terrain configurations. 
Thus, the modeling method must be designed to make good 



Continuous Slope Terrain Profile 




B. Hidden Obstacle Terrain Profile 

FIGURE 10. Possible Terrain Configurations Given 
the Same Two Discrete Sensor Measurements 


FIGURE 11. Illustration of Hidden Passable Terrain 
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"guesses" on the terrain configuration, and, due to the nature 
of the mission, these guesses should be on the conservative 
side. 

A typical situation encountered in the simulation work 
was the occurrence of passable terrain which was hidden from 
the view of the vehicle by a terrain configuration near the 
vehicle's- location. An example of such terrain is shown in 
Figure 11. The terrain presents no actual obstacle but the 
vehicle could not "see" it from the vehicle location shown. If 
the hidden area was large enough so that ^ T^, then 

the modeling method would designate the area as an obstacle 
region. A possible means to overcome this situation would be 
to form a terrain model with two types of unsatisfactory ter- 
rain, namely, confirmed obstacles, obtained by using the slope 
index, vertical index, etc., and potential obstacles, obtained 
by using the discontinuity index. Then, if the vehicle could 
not find a suitable path from its original vantage point, it 
might move in small increments towards a potential obstacle in 
the hope that it would be able to get to a better vantage point. 

Even though the modeling techniques presented in this 
report were designed for long range terrain modeling, it would 
seem that they could be adapted for short range terrain model- 
ing. If it is desired to extend these techniques for this 
case, certain factors should be considered: 

a. The index test numbers should be chosen specifically 
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for the range of terrain to be modeled. 

b. Additional indices may be required to account for 
obstacles peculiar to short range modeling, e.g., 
size and spacing of vertical obstacles. 

c. Small angle approximations that had been used to 
simplify the numerical calculations may not be valid 
for short range modeling, in which case additional 
computing requirements would be necessary to evaluate 
the trigonometric functions. 


VI. CONCLUSION 

The problem considered in this report is the mathematical 
modeling of major terrain features, using measurements acquired 
at a single terrain location, to be used in a self-contained 
surface vehicle guidance system. The terrain modeling problem 
was posed for long-range vehicle guidance, and was conceived to 
consist of three function parts: terrain data acquisition, ter- 
rain data processing, and terrain model generation. 

In Section II, terrain sensor measurements which must pro- 
vide the inputs to any long range terrain modeling system were 
simulated by using a sample terrain, Basic properties of the 
terrain sensor measurements were deduced from the empirical 
data compiled from the sample terrain. 

In Section III, specif ie methods of processing sensor 
measurements were developed in the form of terrain indices. 


The terrain indices would enable the vehicle to detect the 
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presence of an obstacle to travel by processing the acquired 
sensor measurements. A number off indices were developed to 
provide flexibility in designing a terrain modeling system. 

The indices vary in the complexity of the calculations that 
must be performed and the accuracy that can be expected in 
detecting an obstacle. 

In Section IV, complete terrain modeling methods were 
developed using various combinations of the terrain indices 
described in the previous section. These methods were compared 
with respect to their terrain modeling effectiveness and their 
associated computer and sensor requirements. 

This investigation has provided three contributions in the 
area of long range vehicle guidance for an exploratory surface 
vehicle. First, it has developed methods to perform long range 
terrain modeling with only modest computing requirements. 
Second, it has developed general terrain sensing and computer 
requirements necessary to perform long range terrain modeling 
for a selected sample terrain. Third, it has provided general 
techniques to formulate terrain modeling methods for varied 
terrain situations. 

Future work in terrain modeling could include: 1). the 
testing of the modeling methods on other representative ter- 
rains to optimize the performance of the methods, 2). extend- 
ing the techniques developed to the problem of short range 
terrain modeling for surface vehicle guidance. 
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APPENDIX 

Derivation of the Slope Formulas Required 
for the Slope and Gradient Indices 


1. In-Path Slope 

The In-Path Slope (Sj) of a terrain segment is defined as 
the linear approximation to the slope of a terrain segment in 
the direction of a constant azimuth terrain profile, and it is 
expressed as the tangent of the corresponding slope angle, as 
indicated in Figure 12. It can be approximated by dividing the 
change in altitude by the corresponding change in the horizontal 
distance, 

S T - (Al) 

1 A D 


If the altitude Z is described as a function of position in 
polar coordinates, Z = Z (D, 0 ) with the vehicle at the origin, 
then Sj can be interpreted as the numerical approximation to 
the partial derivitive of Z with respect to D, 


Z A 2 

>D AD 


(A2) 


This approximation becomes exact in the limit as A D approaches 
zero. 


The in- path slope can be written in terms of two succes- 
sive sensor measurements acquired along a fixed azimuth 
direction. The quantities R-^ and R 2 denote the corresponding 
range measurements for the elevation angles f3 and 
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( /®2 = /®l + ^Z 3 ) , Sj can be rewritten in terms of the two 


sensor measurements as 


S = R2 sin/g 2 * R 1 sin /^x 

1 R2 cos p 2 - % cosfS-^ 


(A3) 


Since only small values of elevation angle will be used (see 
Section IIC) , the following approximations are made: sin yS ^ ft 

and cos ^3 — 1. Inserting these approximations into equation (A3) , 


Sj becomes 


= r 2 /3 2 - Rj/i 


r 2 


- Ri 


<A4) 


which is in a form that may be easily implemented on a small 
computer . 

2. Cross-Path Slope 

The Cross-Path Slope (S c ) of a terrain segment is defined 
as the linear approximation to the slope of a terrain segment 
in a direction perpendicular to a constant azimuth terrain 
profile, and it is expressed as the tangent of the correspond- 
ing slope angle as indicated in Figure 13. It can be approxi- 
mated by dividing the change in altitude by the corresponding 
change in the horizontal distance (arc length), 


= A 2 

DA O 


(A5) 


If the altitude Z is described as a function of position in 
polar coordinates, Z = Z (D,©) with the vehicle at the origin, 
then S c may be interpreted as the numerical approximation to the 
partial derivative of Z with respect to 0 divided by D, 



FIGURE 13. Illustration of the S r Calculation 
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% Z ^ A Z 
D 2)0 DAO 


(A6) 


where the approximation becomes exact in the limit as A 9 
approaches zero. 

The formula for S c may be written as 


= (z 2 - Z*) 

D 2 a & 


(A7) 


where Z 2 and Z* are discrete points on the terrain each at a 
distance D 2 from the vehicle. If the point Z 2 can be measured 
directly from sensor data, i.e., Z 2 = R 2 sin ^ 2 s at a dis- 
tance D 2 , i.e., D 2 = R 2 cos ^ then, the altitude Z* cannot 
be measured directly from sensor measurements since only ® 
and ^ can be specified by the vehicle's terrain sensor (see 
section IIB) . The altitude Z* can be indirectly measured by 
first computing the in-path slope S-j- from two successive 
sensor measurements Rp j (3 p and R' 2 j|® 2 acc l u i re d along the 

azimuth line O - ©. - A ,® . S' will be an in-path slope 

J ““ 

approximation in the vicinity of Z*, since £3 2 is held fixed 
to acquire R 2 and R' 2 and j3 ^ equals ^3 2 - The slope 

S'j can be computed from these sensor measurements using 
equation (A3) so that 

S' = R*2sin 2 ~ R *1 s ^ n 

1 R'2 cos^3 2 - R'^ cos ^ 1 


(A8) 



which may also be written in terms of Z* as 

S ' j = R *2 sin Pi " Z * (A9) 

R 2 cos ^ 2 ” R 2 cos 2 

This equation can be solved for Z* in terms of sensor measure- 
ments. If the indicated substitutions for Z 2 , Z* and D 2 are 
inserted into (A7) and the small angle approximations for cos p> 
and sin p> are again used, S c can be written as 

s c - (£2 - S ' X >< R 2 - R V (A10) 

R 2 A9 
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When either an in-path or cross-path slope calculation is 
made, the sensor measurements needed for the calculation must 
be acquired from the same terrain feature in order to obtain 
valid slope approximations. A processing system can estimate 
the applicability of a slope calculation for a given set of 
sensor measurements by testing the size of the range differ- 
encing factor (e.g., R 2 - R^) that is required in the slope 


calculation. 
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